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PREFACE 


This  study  was  conducted  as  part  of  the  Independent  Research  Program  of  the 
Institute  for  Defense  Analyses,  under  which  significant  issues  of  general  interest  to  the 
defense  research  community  are  investigated. 
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SCOPE 


The  Monte  Carlo  Layered  Defense  Model  consists  of  six  major  parts.  They 
are  summarized  below. 

1.  Set  up  target  data  base,  including  the  value  of  each  target  and  the  number 
of  terminal  interceptors  at  each  target. 

2.  Set  up  attack  allocation.  Assign  the  RVs  of  all  of  the  attacking  missiles 
to  targets.  Equip  each  missile  with  heavy  and  light  decoys  if  desired. 

3.  Perform  boost-phase  defense,  using  one  of  three  possible  BM/C3 
weapon  assignment  schemes: 

a.  Random  Assignment:  Each  boost-phase  defender  acts  independ¬ 
ently  of  the  others.  The  defender  picks  a  missile  target  and  shoots  at 
it  Thus  each  missile  may  receive  from  none  to  many  shots.  This 
may  be  thought  of  as  unordered  fire,  fully  decentralized  BM/C3,  or 
informationless  fire. 

b.  Efficient  Assignment:  Boost-phase  defenders  are  allocated  as 
uniformly  as  possible  across  their  missile  targets.  For  instance,  if 
there  are  1500  missiles  and  1000  boost-phase  interceptors,  1000  of 
the  1500  missiles,  selected  at  random,  receive  one  shot  each.  If 
there  are  1500  missiles  and  2000  boost-phase  interceptors,  1000 
missiles  selected  at  random  receive  one  shot  each  and  the  other  500 
missiles  receive  two  shots  each.  This  may  be  thought  of  as  ordered 
fire,  centralized  BM/C3,  or  uniform  simultaneous  m-on-n  defense. 

c.  Proportional  Assignment:  Boost-phase  defenders  are  allocated  to 
missiles  in  proportion  to  the  number  of  RVs  on  the  missiles. 

4.  Perform  midcourse  defense,  using  one  of  three  possible  BM/C3  weapon 
assignment  schemes: 

a.  Random  Assignment:  Analogous  to  boost-phase  defense, 
discussed  above. 
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b.  Efficient  Assignment:  Analogous  to  boost-phase  defense, 
discussed  above. 

c.  Preferential  Assignment:  Midcourse  defenders  are  assigned  to  RVs 
to  shape  the  attack  so  that  the  RVs  reaching  targets  chosen  to  be 
defended  are  likely  to  be  killed  by  the  terminal  interceptors.  To 
conserve  the  limited  numbers  of  midcourse  defenders,  the 
algorithm  proceeds  as  follows:  First,  all  targets  receiving  less  than 
a  specified  upper  limit  of  RVs  per  terminal  interceptor  are  defended 
so  that  the  RVs  heading  for  the  targets  are  reduced  below  a 
specified  lower  limit  Second,  the  upper  limit  is  increased  and 
another  set  of  targets  is  defended.  Eventually,  the  lower  limit  is 
reduced  if  there  are  enough  midcourse  interceptors.  The  algorithm 
has  not  been  proven  to  be  optimal,  but  it  is  in  the  spirit  of  one 
known  to  be  optimal  for  adaptive  terminal  preferential  defense. 

(Important  Note:  It  turns  out  that  in  cases  where  the  number  of 
midcourse  interceptors  is  relatively  small  and  damage  is  moderate, 
preferential  assignment  achieves  significant  advantages  over 
efficient  assignment  However,  when  the  number  of  midcourse 
inttictptors  is  increased  and  damage  becomes  very  small, 
preferential  assignment  is  not  superior  to  efficient  assignment). 

Decoys  are  included  in  the  midcoursc  defense  as  follows.  For  any  of 
the  above  three  cases,  when  it  is  decided  to  attack  an  RV  with  an 
interceptor,  the  RV  is  examined  in  the  context  of  all  of  the  objects  from 
its  missile.  The  capability  of  the  defense  to  discriminate  RVs  from 
decoys  is  represented  by  input  probabilities  of  (1)  missed  detection 
against  RVs,  (2)  false  alarm  against  heavy  decoys,  and  (3)  false  alarm 
against  light  decoys.  A  random  process  is  used  to  cither  (1)  attack  the 
RV,  (2)  fail  to  attack  the  RV  because  of  missing  a  detection,  or  (3) 
waste  an  interceptor  on  a  decoy. 

Perform  terminal  defense,  using  any  of  four  possible  BM/C3  weapon 


assignment  schemes: 


a.  Pre  allocated  Fixed-Salvo  Assignment:  This  is  an  implementation  of 
the  Firing  process  associated  with  the  Prim-Read  defense.  This  logic 
is  consistent  with  the  assumption  that  at  the  individual  targets  the 
attacks  are  sequential  and  of  unknown  size.  The  defender  decides  in 
advance  to  fire  a  volley  of  one  or  more  interceptors  at  each  incoming 
RV.  If  one  o*  the  interceptors  in  the  volley  kills  the  RV,  it  does  not 
penetrate.  The  number  of  interceptors  scheduled  to  be  fired  at  all  of 
the  incoming  RVs  adds  up  to  the  total  number  defending  die  target. 
For  instance,  if  there  are  10  interceptors  at  a  target,  preallocated  to 
fire  at  RVs  number  1  through  6,  such  a  firing  pattern  might  be 

3, 2, 2,1, 1,1  with  no  interceptors  fired  at  RV  number  7.  Thus  RV 
number  7  will  surely  kill  the  defended  target,  and  any  of  the 
engagements  against  RVs  1  through  6  might  also  result  in  the 
defended  target  being  killed. 

b.  Efficient  Assignment:  Analogous  to  boost-phase  defense  and  mid¬ 
course  defense,  discussed  above.  Consistent  with  simultaneous 
attacks  of  known  size  at  the  individual  targets. 

c.  Limited  Shoot-Look-Shoot:  The  defender  has  a  certain  number  of 
interceptors  at  the  target  At  the  first  RV  he  fires  interceptors  one- 
by-one  until  either  the  RV  is  killed  or  the  upper  bound  on 
interceptors  able  to  be  fired  is  reached.  For  instance,  there  might  be 
up  to  three  shots  at  each  RV.  If  all  three  fail  to  kill  the  RV  there  is  a 
penetrator.  This  doctrine  conserves  interceptors,  unlike  ihe  two 
doc  rimes  above,  for  there  arc  no  interceptors  "wasted"  in  salvos. 

d.  Unlimited  Shoot-Look-Shoot:  Same  as  limited  shoot-look-shoot 
above,  but  no  upper  bound  on  number  of  shots  at  each  RV.  This 
doctrine  produces  the  greatest  possible  effectiveness  of  a  given 
number  of  terminal  interceptors  at  a  defended  target. 

6.  Assess  damage  to  targets.  For  each  target,  determine  if  there  are 
one  or  more  penetrating  RVs. 

There  are  two  significant  features  of  the  model  available  to  extend  and 
elaborate  on  the  above,  as  follows: 
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1.  It  is  possible  to  suppress  the  boost-phase  defense,  either  by  killing  a 
certain  number  of  defenders  or  by  permitting  a  certain  r, umber  of 
missiles  not  to  be  engaged  by  the  boost-phase  defense. 

2.  It  is  possible  to  perform  hybrid  combinations  of  BM/C3  rules  within 
boost-phase  defense  and  midcourse  defense.  Specifically,  a  certain 
portion  of  the  missiles  in  boost-phase  defense,  or  RVs  in  midcourse 
defense,  can  be  engaged  efficiently  by  a  portion  of  the  defense  while  the 
remainder  are  engaged  randomly  by  a  portion  of  the  defense.  This 
permits  exploration  of  a  partially  centralized  engagement  followed  by  a 
partially  decentralized  engagement  of  the  survivors. 


KEY  ASPECTS 


Three  key  aspects  of  the  qualitative  behavior  of  the  model  are  summarized  below. 

1.  Random  Versus  Efficient  Weapon  Allocation  in  Boost-Phase  and 
Midcourse  Defense 

Experience  with  the  model  has  shown  that  there  are  very  significant 
effects  associated  with  random  versus  efficient  assignment  of  boost- 
phase  defenders.  Completely  decentralized  assignment  of  boost-phase 
defenders  leads  to  many  missiles  not  being  attacked,  and  thus 
substantially  more  RVs  being  presented  to  the  midcourse  defense. 

Experience  with  the  model  has  also  shown  that  there  are  very 
significant  effects  associated  with  random  versus  efficient  assignment  of 
midcourse  defenders.  If  midcourse  defenders  can  be  efficiently 
assigned,  at  a  certain  point  every  single  RV  is  confronted  by  a  one-on- 
one  engagement,  or  a  two-on-one  engagement,  and  thus  there  tend  to  be 
few  cases  in  which  the  terminal  defenses  are  exhausted. 

The  amount  of  resources  required  to  achieve  near- zero  damage  is 
extremely  sensitive  to  the  BM/C3  logic.  Random  assignment  of  boost- 
phase  and  midcourse  defenders  leads  to  twice  as  many  boost-phase 
defenders  and  midcourse  defenders  being  required  to  reduce  damage  to 
near- zero  levels  than  does  efficient  assignment,  in  a  test  case. 

To  achieve  efficient  defense  in  boost-phase  and  midcourse  requires 
that  the  defender's  weapons  be  appLied  to  missiles  and  to  RVs 
sparingly.  First  consider  boost-phase  defense.  Assume  that  each 
missile  is  to  be  engaged  once.  Assume  that  there  is  no  shoot-look- 
shoot.  If  the  attack  is  simultaneous,  the  centralized  system  must  allocate 
all  of  the  shooters  to  all  of  the  targets  uniformly.  This  must  be  done 
within  a  period  of  a  minute  or  so.  If  the  attack  is  sequential,  say  over 
five  or  ten  minutes,  the  centralized  system  must  allocate  the  shooters  to 
the  subset  of  targets  presented  in  accordance  with  a  preallocated 
rationing  scheme.  If  one-half  of  the  missiles  have  been  launched,  one 
half  of  the  defenders  should  be  utilized  (or  something  similar  to  this  ). 


Partitioning  the  boost-phase  defenders  into  subsets  should  not 
preclude  efficient,  as  opposed  to  random,  defense  if  the  subsets  are 
large  enough.  Each  subset  would  be  responsible  for  a  subset  of  targets. 
Difficulties  in  efficient  assignment  would  occur  if  the  shooters  and 
targets  could  not  be  related  one-to-one.  Since  the  shooters  orbit  out  of 
range  of  targets,  time-phased  assurance  that  all  missile  targets  receive 
about  the  same  number  of  shots  may  be  difficult,  but  should  be  stressed 
as  a  design  parameter  of  the  boost-phase  defense. 

Now  consider  midcoursc  defense,  which  takes  place  over  a 
relatively  long  period  of  time.  By  monitoring  each  shooter  and  RV 
target,  it  should  be  possible  to  assure  that  each  RV  receives  about  the 
same  number  of  shots.  With  respect  to  time  considerations,  efficient 
allocation  seems  easier  in  midcourse  than  in  boost-phase.  Simultane. 
attack  versus  ragged  attack  does  not  seem  to  have  significant  timing 
effects  in  midcourse  defense  from  an  assignment  point  of  view. 

The  principal  problem  in  midcourse  defense  seems  to  be  the  effect  of 
large  numbers  of  decoys.  If  the  K-factor,  which  is  basically  the  number 
of  standard  deviations  of  the  measuring  process  separating  RVs  from 
decoys,  is  large,  then  few  RVs  will  escape  detection  and  few 
interceptors  will  be  wasted  on  decoys.  If  the  K-factor  is  small,  or  if 
there  are  enough  light  decoys  with  even  a  small  probability  of  being 
engaged  (say  0.1),  then  the  midcourse  interceptors  may  be  exhausted  by 
the  decoys,  and  the  efficient  allocation  of  defenders  to  RVs  will  fail. 

-■  Terminal  Defense 

The  function  of  the  terminal  defense  is  to  cope  with  the  RVs  which 
get  through  to  attack  the  target  The  number  of  RVs  may  vary  widely, 
since  if  there  are  many  missiles  shooting  at  the  target  there  may  be  any 
number  from  zero  to  a  large  number  of  RVs  heading  for  the  target  after 
the  boost-phase  defense.  The  first  two  types  of  midcourse  defense  will 
not  eliminate  this  phenomenon.  The  third  type  will  solve  it  for  some 
targets,  but  othei  targets  will  be  heavily  attacked  since  they  are  "written 
off'  by  the  rrudeourse  defense. 
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Performance  of  the  terminal  defense  in  situations  where  the  number  of 
defenders  is  just  a  bit  larger  than  the  number  of  attacking  RVs  provides 
the  key  difference  in  effectiveness.  For  instance,  if  there  are  seven 
defenders  with  P^.9,  and  six  attackers,  a  preallocated  defense  with 

Firing  doctrine  2, 2, 1,1, 1,0  would  yield  no  chance  of  survival,  while 
shoot-look-shoot  would  yield  a  50  percent  chance  of  survival.  The 
Monte  Carlo  process  results  in  quite  different  engagements  at  targets,  of 
which  a  significant  number  are  in  the  regions  where  firing  doctrines 
make  a  big  difference. 

3.  Expected  Value  Versus  Monte  Carlo 

Another  key  aspect  of  the  overall  performance  of  the  model  is  in  the 
difference  between  an  expected  value  model  and  a  Monte  Carlo  model.  A 
typical  expected  value  model  makes  an  assessment  of  the  expected 
number  of  missiles  killed  by  the  boost-phase  defense  and  then  counts 
the  expected  number  of  RVs  left  It  assumes  that  these  RVs  are  targeted 
proportional  to  the  targeting  of  the  missiles  from  whence  they  came. 
Then  the  expected  value  model  makes  an  assessment  of  the  number  of 
RVs  killed  by  the  midcourse  defense.  The  RVs  surviving  are  again 
assumed  to  be  targeted  proportional  to  the  original  targeting  of  the 
missiles.  These  RVs  are  engaged  by  the  terminal  defenders.  Typically, 
the  terminal  defenders  significantly  outnumber  the  RVs  at  all  of  the 
targets  because  the  RVs  are  uniformly  allocated  to  the  targets,  and  thus 
only  one  or  two  RVs  per  target  are  confronted  by  at  least  one  interceptor 
each.  The  Monte  Carlo  model,  on  the  other  hand,  results  in  a  significant 
variations  in  RVs  appearing  at  the  attacked  targets.  This  one  property  is 
lhe_main  rcason_for  the  necessity  for  the  Monte  Carlo  model.  Some 
experimental  results  have  shown  that,  when  the  of  boost-phase, 

midcourse  and  terminal  defense  is  .9,  the  expected  percent  of  the  data 
base  destroyed  in  the  Monte  Carlo  Layered  Defense  Model  is  5. 1 
percent  for  preallocated  terminal  defense  and  2.1  for  efficient  terminal 
defense,  whereas  the  expected  percent  of  the  data  base  destroyed  by  an 
expected  value  calculation  is  1.0  percent  and  .3  percent,  respectively. 


For  near-pcrfect  defense,  or  the  assured  survival  mission  of  SDI,  these 
differences  are  crucial. 
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SUBROUTINES 


a.  The  main  indexes  are  II,  12,  and  13,  where 

(1)  II  indexes  on  the  number  of  boost-phase  interceptors 

(2)  12  indexes  on  the  number  of  midcourse  interceptors 

(3)  13  indexes  on  the  number  of  combinations  of  boost-phase 
and  midcourse  BM/C3  options 

b.  For  each  value  of  13,  three  indicators  an*,  set: 

(1)  131 -indicator  for  one  of  three  boost- phase  defense  options 

(2)  132— indicator  for  one  of  three  midcourse  defense  options 

(3)  133— indicator  for  one  of  four  terminal  defense  options 

The  options  are  discussed  later.  In  addition,  the  indicators  IB  PM 
and  IMCM  are  set: 

(1)  IBPM— indica  tor  for  hybrid  boost-phase  defense,  where  0 
means  no  hybrid  and  1  means  hybrid 

(2)  IMCM— indicator  for  hybrid  midcourse  defense,  where  0  means 
no  hybrid  and  1  means  hybrid 

c.  The  long  print  indicator  LPI  is  set: 

(1)  LPI=0  means  no  long  printout 

(2)  LPI“  1  means  long  printout 

d.  The  random  number  seed  IRS  for  the  entire  simulation  is  set-  The 
seed  is  in  COMMON  and  gives  a  new  value  each  time  a  random 
number  is  generated  anywhere  in  the  model. 

e.  The  value  target  data  base  is  set  up  by  calling  Subroutine 

VALUES(NVT).  Subroutine  VALUES  creates  array  IV (I), 

1=1, NVT  consisting  of  the  values  of  targets  1  through  NVT.  It 
provides  NVT  to  MAIN.  Array  IV  is  in  COMMON. 

f.  The  main  loop  on  the  number  of  replications,  irdexed  by  IMAIN 
from  1  through  NUMB,  is  begun  by  the  following  statement: 
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DO  8000  MAIN  -  1,  NUMB 

g.  The  terminal  defense  is  set  up  by  calling  Subroutine  TDLNV  (NVT). 
Subroutine  TDINV  creates  array  ITD  (I),  1=1  NVT  consisting  of  the 
number  of  terminal  interceptors  at  targets  1  through  NVT.  Array 
ITD  is  in  COMMON. 

h.  The  attack  allocation  is  set  up  by  calling  Subroutine  ATT  ALL 

(NMIS,  MNWT).  Subroutine  ATT  ALL  creates  arrays  as  follows: 


IAT  (U), 

I=i,NMIS, 

J=1MNWT 

LAW  (I,  J), 

1=1  NMIS, 

J-1.MNWT 

IAD(U,1) 

I=J,NMIS, 

1-1  MNWT 

IAD(1,J,2) 

1=1, NMIS, 

J=1,MNWT 

NMIS  is  the  total  number  of  missiles  and  MNWT  is  the  maximum 
number  of  targets  which  each  missile  attacks.  For  instance,  NMIS 
would  be  1500  and  MNWT  would  be  2  when  there  are  1500 
missiles,  each  attacking  no  more  than  two  different  targets. 

IAT(IJ)  contains  the  target  or  targets  to  which  the  RVs  of  missile 
I  are  assigned  For  instance,  if  missile  10  is  assigned  to  targets  100 
and  200,  then  IAT(  10,1)=  100  and  IAT(10,2)=200. 

IAW(IJ)  contains  the  number  of  RVs  assigned  from  missile  I  to 
the  targets  specified  by  IAT(I,J).  In  the  above  example,  if  missile 
10  has  8  RVs,  and  if  four  each  are  assigned  to  targets  100  and  200, 
then  LAW  (10,1)  =  4  and  IAW(10,2)  =  4. 

LAD  (L  J,  1)  and  IAD  (I,  J,  2)  contain  the  number  of  decoys  from 
missile  I  to  the  targets  specified  by  LAT  (I,  J).  In  the  above 
example,  if  missile  10  has  4  heavy  decoys  and  50  light  decoys  split 
evenly  between  the  two  targets  then  IAD  (10,  1,  1)  =  2,  IAD  (10,  1, 
2)  =  25,  IAD  (10,  2,  1)  =  2  and  IAD  (10,  2,  2)  -  25. 


Arrays  I  AT,  EAW  and  LAD  are  in  COMMON.  The  values  of  NMIS 
and  MNWT  are  specified  in  Subroutine  ATT ALL  and  returned  to 
MAIN. 

The  boost-phase  defense  is  activated  against  the  missiles  by 
calling  Subroutine  BPD(NSHT,  NTAR,  MNWT,  PK,  END), 
where: 

NSHT  is  the  number  of  boost-phase  defenders 

NTAR  is  the  number  of  missile  targets  (typically  set 
equal  to  NMIS) 

MNWT  is  the  maximum  number  of  value  targets 
attacked  by  each  missile. 

PK  is  the  boost-phase  defense  single-shot  probability  of  kill. 

END  is  the  defense  option,  as  fellows: 

(1)  END  =  1  Random  allocation  of  boost-phase  defenders.  Each 
boost-phase  defender  selects  and  fires  at  missile  independently  of  the 
other  defenders. 

(2)  END  -  2  Efficient  allocation  of  boost-phase  defenders.  Boost- 
phase  defenders  are  allocated  uniformly  over  missiles. 

(3)  END  -  3  Proportional  allocation  of  boost-phase  defenders. 
Boost-phase  defenders  are  allocated  to  missiles  in  proportion  to  the 
number  of  RVs  on  the  missiles. 

The  indicator  for  the  hybrid  boost-phase  defense  EBPM  is  queried 
before  Subroutine  9 PD  is  called: 

(1)  If  IBPM  =  0,  then  the  number  of  shots  equals  NSHT,  as 
above. 

(2)  If  EBPM  =  1,  then  the  boost  phase  engagement  proceeds  in  two 
steps.  In  the  first  step  a  certain  number  of  shoots  is  assigned  to 
efficient  defense  (IND  =  2).  In  the  second  step  the  remainder  is 
assigned  to  random  defense  (END  =  1). 


j.  The  midcourse  defense  is  activated  against  the  RVs  by  calling 
Subroutine  MOD  (NVT,  NSHT,  NTAR,  MNWT,  PK,  IND),  where: 

NVT  is  the  number  of  value  targets  being  protected 

NSHT  is  tine  number  of  midcourse  defenders 

NTAR  is  the  number  of  missiles  from  which  the  RVs  are 
being  released  (typically  set  equal  to  NMIS) 

MNWT  is  the  maximum  number  of  value  targets  attacked  by  each 
missile 

PK  is  the  midcourse  defense  single-shot  probability  of  kill 
IND  is  the  defense  option  as  follows: 

(1)  IND  =  1  Random  allocation  of  midcourse  defenders.  Each 
midcourse  defender  selects  and  fires  at  RVs  independently  of  the 
other  defenders. 

(2)  IND  =  2  Efficient  allocation  of  midcourse  defenders.  Midcourse 
defenders  are  allocated  uniformly  over  RVs. 

(3)  IND  =  3  Preferential  allocation  of  midcourse  defenders. 
Assuming  knowledge  of  the  complete  attack,  including  the 
destination  of  all  RVs,  midcourse  interceptors  are  allocated  to  shape 
the  attack  such  that  the  RVs  which  remain  after  the  midcourse  defense 
can  be  destroyed  by  the  terminal  defense,  the  overall  goal  being  to 
preserve  as  much  surviving  value  as  possible. 

The  indicator  for  the  hybrid  midcourse  defense  IMCM  is  queried 
before  Subroutine  MCD  is  called: 

(1)  If  IMCM  =  0,  then  the  number  of  shooters  =NSHT,  as  above. 

(2)  If  IMCM  *  1,  then  the  midcourse  engagement  proceeds  in  two 
steps.  In  the  first  step  a  certain  number  of  shooters  is  assigned  to 
efficient  defense  (IND  =  2).  In  the  second  step  the  remainder  is 
assigned  to  random  defense  (IND  =1). 

k.  The  terminal  defense  is  activated  against  RVs  arriving  at  each  target 

by  calling  Subroutine  TD  (NVT,  NMIS,  MNWT,  PK,  IND),  where 
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N  VT  is  the  number  of  value  targets  being  protected 

NMIS  is  the  number  of  missiles  in  the  attack 

MNWT  is  the  maximum  number  of  value  targets  attacked  by  each 
missile 

PK  is  the  terminal  defense  single-shot  probability  of  kill 
IND  is  the  defense  option,  as  follows: 

(1)  IND  =  1  Freallocated  fixed-salvo  defense.  For  each  target,  each 
arriving  RV,  sequentially,  receives  a  salvo  of  a  pre-specified  number 
of  interceptors. 

(2)  DSHD  =  2  Efficient  allocation  of  terminal  interceptors.  For  each 
target,  interceptors  are  allocated  as  uniformly  as  possible  over  all 
attacking  RVs. 

(3)  IND  =  3  Limited  shoot-look-shoot.  For  each  target,  RVs  are 
engaged  one-by-one.  The  First  RV  is  shot  at  by  one  interceptor,  then 
another,  until  either  it  is  killed  or  the  upper  limit  per  RV  is  reached. 
The  second  RV  is  engaged  similarly.  The  process  terminates  when 
there  are  no  more  RVs  or  no  more  interceptors. 

(4)  IND  -  4  Unlimited  shoot-look-shoot.  Same  as  (3)  above,  but 
no  limitation  on  interceptor  engagements  per  RV. 

l.  After  the  terminal  defense  there  are  NW(I),  I  =  1,NVT  succcessful 
penetrators  at  each  target  Where  NW(I)  =  0  the  target  is  not 
destroyed.  The  terminal  defense  subroutine  TD  creates  the  relative 
frequency  array  of  successful  penetrators  NWRF  (I). 

m.  Value  surviving  after  the  attack  is  computed  by  Subroutine 
VALSURV  (NVT,  IVSS),  which  creates  the  value  sum  surviving 
associated  with  those  targets  where  NW  (I)  =  0, 1  =  1,  NVT. 

n.  The  average  and  standard  deviation  of  value  damaged  for  the  trials 
made  thus  far  (up  to  the  current  value  of  IMAIN)  is  computed  and 
printed. 


o.  This  is  the  end  of  the  main  loop  numbered  8000  on  the  number  of 
iterations  to  be  performed. 

p.  Tills  is  the  end  of  the  main  loop  numbered  9999  on  the  specif  ication 
of  the  case  to  be  simulated,  for  a  particular  number  of  boost-phase 
defenders  and  midcourse  defenders,  and  for  a  defense  logic  for  the 
three  layers. 

2.  Subroutine  VALUES  (NVT) 

Subroutine  VALUES  (NVT)  creates  the  array  of  target  values  IV, 
where  IV(I)  is  the  value  of  target  I,  and  returns  to  MAIN  the  quantity 
NVT,  the  total  number  of  targets.  The  content  of  the  subroutine  is  at  the 
option  of  the  user. 

For  instance,  if  there  were  1000  targets  with  values  1000  through  1, 
the  contents  of  the  subroutine  would  be  as  follows: 

NVT  =  1000 

DO  10  I  =  1,  1000 

IT  =  1001  -  I 

10  IV(I)  IT 

Subroutine  VALUES  (NVT)  of  the  example,  contained  in  Appendix 
B,  generates  a  set  of  values  for  2533  targets. 


3.  Subroutmg-JD.INy. (NVD 


Subroutine  TDINV  (NVT),  Terminal  Defense  Inventory,  creates  the 
array  of  terminal  inventories  of  interceptors  ITD,  where  TIT)  (I)  is  the 
number  of  terminal  interceptors  at  target  I.  The  content  of  the 
subroutine  is  at  the  option  of  the  user. 


For  instance,  if  in  subroutine  VALUES  there  were  created  1000 
targets  with  values  1000  through  1,  and  if  they  were  to  be  defended  by 
one  interceptor  per  ten  units  of  value,  truncated  to  the  nearest  integer, 
die  contents  of  the  subroutine  would  be  as  follows: 
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DO  10  1=1,  1000 
10  ITD  (I)  =  IV(I)  /  io 


Subroutine  TDINY  (NVT)  of  the  example,  contained  in  Appendix 
B,  deploys  terminal  interceptors  for  2533  targets  equal  to  one-fourth  of 
their  value,  truncated  to  the  nearest  integer. 


4.  Subroutine  ATT ALL  (NMIS.  MNWTl 

Subroutine  ATrALL(NMIS,  MNWT),  Attack  Allocation,  contains 
the  allocation  of  attackers  to  targets. 

It  first  sets  the  number  of  missiles  NMIS  and  the  maximum  number 
of  targets  attacked  by  each  missile  MNWT  (equivalently,  maximum 
number  of  targets  cross-targeted  by  any  attacking  missile). 

The  subroutine  creates  arrays  1AT  (I J),  LAW  (IJ),  IAD  (I,J,K) 
where  I  =  1,  NMIS,  J  =  1,  MNWT  and  K  =  1,2. 


Consider  a  simple  example  of  1000  missiles  attacking  targets  1 
through  1000  with  5  RVs  on  each  missile,  2  heavy  decoys  on  each 
missile  and  10  light  decoys  on  each  missile.  Each  missile  attacks  only 
one  target.  The  contents  of  the  subroutine  would  be  as  follows: 
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NMIS  =  1000 
MNWT = 1 . 

DO  10  1=  1,  1000 
IAT  (I,  1)  -  I 
IAW  (I,  1)  =  5 
IAD  (I,  1,  1)  =  2 
IAD  (I,  1,2)=  10 


Subroutine  ATT  ALL,  (NMIS,  MNWT)  of  the  example,  contained  in 
Appendix  B,  sets  up  an  attack  on  500  targets.  It  allocates  500  ICBMs 
with  5  RVs  each  to  the  500  targets,  then  500  ICBMs  with  5  RVs  each  to 
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the  same  500  targets,  then  500  SLBMs  with  10  RVs  each  to  the  same 
500  targets.  In  the  option  where  decoys  are  included,  the  ICBMs  have 
3  RVs,  4  heavy  decoys  and  30  light  decoys  while  the  SLBMs  have  6 
RVs,  8  heavy  decoys  and  30  light  decoys. 

5.  Subroutine  BPD  (NSHT.  NTAR,  MNWT,  PK.  iND 

Subroutine  BPD  (NSHT,  NTAR,  MNWT,  PK,  IND),  Boost  Phase 
Defense,  receives  from  MAIN  the  number  of  boost-phase  shooters 
NSHT,  number  of  missile  targets  NTAR,  maximum  number  of  targets 
attacked  by  any  missile  MNWT,  probability  of  kill  PK  and  indicator  of 
defense  option  END. 

Its  function  is  to  take  the  array  of  missiles  accounted  for  in  IAT(I,  J) 
and  to  kill  certain  of  them.  For  those  killed,  the  arrays  in  COMMON 
IAT  (I,  J),  LAW  (I,  J),  IAD  (I,  J,  1)  and  IAD  (I,  J,  2)  are  set  to  zero. 

a.  Boost-Phase  Defense  Suppression 

The  firs  t  step  is  to  decide  whether  or  not  to  suppress  the  boost- 
phase  defense,  and,  if  so,  how  it  is  to  be  suppressed.  The  indicator 
IBPDS  is  consulted,  with  options  as  follows: 

IBPDS  =0  no  defense  suppression 

IBPDS  =  1  reduce  number  of  shooters  by  NKLD 

IBPDS  =  2  reduce  number  of  shooters  by  NSHT  *FKLD 

IBPDS  =  3  do  not  engage  number  of  missiles  NTARS 

EBPDS  =  4  do  not  engage  fraction  of  missile'-  FTARS 

In  the  beginning  of  the  subroutine  the  user  sets  values  for  IBPDS, 
NKLD,  FKID,  NTARS,  FTARS.  These  are  used  as  necessary', 
depending  on  the  option  indicated  by  IBPDS. 

If  IBPDS  =  0,  1  or  2,  then  IBPDST  =  0.  For  these  three  cases  the 
number  of  boost-phase  shooters  is  reduced  by  0,  NKLD  or 
NSHT*FKLD. 
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If  IBPDS  =  3  or  4,  then  EBPDST  =  1.  The  fraction  of  missiles  not 
to  be  engaged  is  set  to  FTARS. 


b.  Random  Defense 

Consider  boost  phase  defenders  1  through  NSHT.  For  each 
defender,  select  a  random  number  A.  If  IBPDST  -  0,  the  missile 
corresponding  to  that  number  will  be  engaged.  If  IBPDST  =  1 ,  the 
missile  corresponding  to  that  number  will  not  be  engaged  if  A  £ 
FTARS. 

If  the  missile  is  to  be  engaged,  draw  a  random  number  to  see  if  it 
is  killed.  If  it  is  killed,  set  IAT,  IAW  and  LAD  arrays  equal  0  for 
that  missile. 

c.  Efficient  Defense 

(1)  Consider  the  case  where  boost-phase  shooters  are  fewer  than 
missile  targets  (NSHT  <  NTAR).  Set  INDTA  (I),  1=1,  NTAR  - 
0  to  keep  track  of  which  targets  have  not  been  engaged.  Consider 
boost-phase  defenders  1  through  NSHT.  If  IBPDST  =  1,  draw  a 
random  number  A  and  do  not  engage  missile  if  A  ^  FTARS.  If  a 
missile  is  to  be  engaged,  draw  a  random  number  A  to  see  which 
missile  is  to  be  engaged.  Call  this  IMIS.  If  INDTA  (IMIS)  =  1 
skip  this  missile  for  it  has  already  been  engaged.  Instead,  draw 
another  random  number  to  attempt  to  find  a  target.  If  INDTA 
(IMIS)  =  0  then  set  INDTA  (IMIS)  =  1,  draw  a  random  number 
and  determine  if  the  missile  is  killed.  If  it  is  killed,  set  IAT,  IAW 
and  IAD  to  zero  for  that  missile. 

(2)  Now  consider  the  case  where  boost-phase  shooters  are  equal 
to  missile  targets  (NSHT=NTAR).  Do  not  bother  to  set  INDTA 
(I),  1=1,  NTAR  to  0,  since  all  targets  will  be  engaged  unless 
defense  is  suppressed. 

For  boost-phase  defenders  1  through  NSHT,  if  IBPDST  =  1, 
draw  a  random  number  A  and  do  not  engage  missile  if  A  <; 

FT  ARS.  If  missile  is  to  be  engaged,  draw  a  random  number  A  to 
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sec  if  it  is  killed.  If  so,  set  IAT,  IAW  and  IAD  to  zero  for  that 
missile. 

(3)  Now  consider  the  case  where  boost-phase  shooters  are  greater 
than  number  of  missile  targets  (NSHT  >  NTAR).  Set  INDTA  (I),  I 
=  1,  NTAR  equal  to  0  to  keep  track  of  targets  which  have  been 
engaged. 

A  fraction  PI  of  the  targets  will  receive  IRP1  shots  and  a  fraction 
(1  -  PI)  of  the  targets  will  receive  IR  shots,  where 

R  .  SHT/TAR 

IR  -  R 

FIR=  IR 

PI  =R  -  FIR 

IRP1  =  IR  +  1 

For  example,  if  SHT  =  1400  and  TAR  =  1000,  then  .4  of  the 
targets  will  receive  2  shots  and  .6  of  the  targets  will  receive  1  shot. 

First,  consider  the  NTP1  targets  which  will  receive  IRP1  shots. 
Check  to  see  if  IBPDST=1.  If  so,  and  if  A  ^  FTARS,  a  missile  is 
not  attacked.  Otherwise,  pick  a  missile  IMIS  at  random.  If 
INDTA(IMIS)=0  it  can  be  attacked  Set  INDTA(I)=1  and  shoot 
IRP1  shots  at  the  missile.  If  it  is  killed  set  IAT,  IAW  and  IAD  to 
zero  for  that  missile. 

Second,  consider  the  remainder  of  the  targets.  If  ENDTA(I)=0  and 
defense  not  suppressed,  shoot  IR  shots  at  the  missile.  If  it  is 
killed  set  IAT,  IAW  and  IAD  to  zero  for  that  missile. 

d.  Defense  in  Proportion  to  Number  of  RYs 

Count  the  total  number  of  RVs  FIAWT.  Compute  the  desired 
shots  per  warhead  SPW. 
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cv. 


(1)  Consider  the  case  where  shooters  are  less  than  targets  (NSHT 
<  NTAR).  Set  INDTA(I)-0, 1=1,  NTAR  to  0  to  keep  tra^k  of 
targets  which  have  been  engaged. 

Pick  a  missile  at  random  and  count  the  number  of  RVs  on  it 
FIAWMT.  The  number  of  shots  to  be  fired  at  that  missile  is 
SP  W*FIA  WMT,  rounded  to  the  next  higher  or  lower  integer  (using 
a  random  number  based  on  the  number  of  shooters  and  targets). 
Shoot  up  to  that  many  shots  at  the  missile  and  if  it  is  killed  set  IAT, 
LAW  and  I  AD  to  zero. 

(2)  Consider  the  case  where  shooters  are  equal  to  or  greater  than 
targets  (NSHT  ^  NTAR).  There  are  two  options  within  this  case. 

If  at  least  one  shot  per  missile  (complete  coverage)  is  not  required 
ICMPLT-0)  then  fire  shots  at  missiles  in  accordance  with  their 
number  of  RVS.  In  this  case  a  missile  may  not  be  engaged  it  it  has 
fewer  RVs  than  the  cutoff. 

Alternatively,  pass  through  all  of  the  missiles,  firing  one  shot  at 
each.  Then,  utilizing  the  remainder,  fire  in  proportion  to  the 
originally  computed  shots  per  warhead. 

6.  Subrouting  MCE  (NVT^SHTuNIARMWf  ,PK>IND) 

Subroutine  MCD  (NVT,  NSHT,  NTAR,  MNWT,  PK,  IND), 
Midcourse  Defense,  receives  from  MAIN  the  number  of  value  targets 
NVT,  number  of  midcourse  shooters  NSHT,  number  of  RV  targets 
NTAR,  maximum  number  of  value  targets  attacked  by  any  missile 
MNWT,  probability  of  kill  PK  and  indicator  of  defense  option  IND. 

The  indicator  EDECOYS  is  set  at  the  beginning  of  MCD.  If 
IDECOYS=0  there  are  no  decoys.  If  IDECOYS=  1  there  are  decoys. 
Associated  with  decoys  are  parameters  as  follows: 

PMD  =  probability  of  missed  detection 

PFA1  -  probability  of  false  alarm  against  decoys  of 
type  1 
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FFA2  -  probability  of  false  alarm  against  decoys  of 
type  2. 

The  array  IAWT  (IJ)  is  created.  It  is  set  equal  to  IAW  (I, J).  The 
array  IAWT  will  contain  the  actual  number  of  RVs  as  the  subroutine 
progresses.  In  making  decisions  some  of  the  options  will  consult  the 
array  IAW  (I,J),  the  RVs  entering  the  subroutine,  while  others  will 
consult  the  array  IAWT  (I,J),  the  current  status  of  the  simulation. 

a.  No  Decoys.  Random  Defense 

For  shots  1  through  NSHT,  select  an  RV  by  drawing  a  random 
number,  picking  a  missile  and  seeing  if  there  were  any  RVs  on  that 
missile  before  it  entered  subroutine:  MCD.  This  is  done  by 
examining  the  array  IAW  (I,  J). 

If  there  is  one  or  more  RV  associated  with  die  missile,  shoot  at 
the  RV.  If  it  is  killed,  reduce  the  actual  number  in  IAWT  (I,  J) 
down  as  far  as  zero. 

In  this  option  RVs  can  be  killed  even  if  already  dead,  since  there 
is  no  coordination. 

b.  No  Decoys,  Efficient  Defense 

First,  consider  shooters  equal  to  or  less  than  targets  (NSHT  <; 

N TAR).  For  each  shot,  select  a  missile.  If  that  missile  actually 
contains  an  RV  as  indicated  by  IAWT  (I,  J),  shoot  at  it,  and  if 
successful,  reduce  the  actual  RVs  on  the  missile  by  1.  If  that 
missile  does  not  contain  an  RV,  go  on  to  consider  another  missile 
by  again  drawing  a  random  number. 

Next,  consider  shooters  greater  than  targets  (NSHT  >  NTAR). 

In  this  case  plan  to  shoot  at  a  fraction  PI  of  the  targets  IRP1  shots 

and  at  a  fraction  (1 -PI)  of  the  targets  IR  shots,  where,  as  in  the 

boost-phase  defense, 


R  =  SHTTAR 


For  each  missile,  determine  whether  or  not  it  has  any  RVs.  II'  it 
does,  count  the  total  IAWTT.  For  each  of  the  IAWTT,  by 
drawing  a  random  number,  decide  if  it  will  receive  IRP1  or  IR 
shots.  Shoot  up  to  the  desired  number  of  shots  at  the  RV,  If  it  is 
killed,  reduce  IAWT  (I,  J)  by  1. 

.  No  Decovs.  Preferential  Defense 

For  this  case  Subroutine  MCD3  (NVT,  NSHT,  NTAR, 
MNWT,  PK)  is  called. 

The  number  of  RVs  assigned  to  each  target  is  computed  as 
NWT  (IMAIN),  IMAIN=1,  NVT.  For  each  target  indexed  by 
IMA  IN,  the  indicator  ENDI  (IMAIN,  JT)=I  indicates  the  number 
of  the  Ith  missile  shooting  at  it  and  the  indicator  INDJ  (IMAIN, 
IT)  =  J  indicates  the  Jth  group  of  RVs  from  the  Ith  missile.  In  the 
example  there  are  three  missiles  shooting  at  each  target,  and  2533 
targets,  so  the  dimensions  are  XNDI  (2533,3)  and  INDJ  (2533,3). 

Next,  there  are  seven  iterations  wherein  the  midcourse  defense 
attempts  to  shape  the  entire  attack  so  that  the  terminal  defenders 
can  cope  with  it.  Upper  bounds  RATSIJ  and  lower  bounds 
RATSL  are  set  for  the  quantity 

RAT  =  N  VVT(INlAlN)/nTD(IMAIN) 

or  ratio  of  attackers  to  terminal  defenders.  The  least-attacked 
targets  are  defended  first  (as  measured  by  an  upper  bound', 
reducing  the  RVs  until  they  ar::  below  a  certain  lower  bound. 
Then  another  set  of  controls  is  established,  and  so  on.  The 
algorithm  is  summarized  as  follows: 


Count  number  of  RVs  assigned  to  each  target  from  all  surviving 
missiles.  At  each  target  let  X=RVs, Terminal  Defenders.  Until 
midcourse  defense  is  exhausted: 

(1)  For  all  targets  where  X  5  1  kill  sufficient  RVs  to  achieve 
X  £  .8 

(2)  For  all  targets  where  X  £  2  kill  sufficient  RVs  to  achieve 
X  5.8 

(3)  For  all  targets  where  X  5  3  kill  sufficient  RVs  to  achieve 
X  5.8 

(4)  For  all  targets  where  X  5  4  kill  sufficient  RVs  to  achieve 
X  5  .8 

(5)  For  all  targets  where  X  5  U*  kill  sufficient  RVs  to  achieve 
X  5  .8 

(6)  For  all  targets  where  X  5 .8  kill  sufficient  RVs  to  achieve 
X  5.6 

(7)  For  all  targets  where  X  5  .6  kill  sufficient  RVs  to  achieve 
X  =  0, 


where  U*  is  the  input  upper  bound  (currently  10). 

As  RVs  are  killed  they  are  subtracted  from  IAWT. 
d.  Decovs.  Random  Defense 

For  shots  1  through  NSHT,  an  RV  is  selected  by  drawing  a 
random  number,  picking  a  missile  and  seeing  if  there  were  any 
RVs  on  that  missile  before  it  entered  subroutine  MCD.  This  is 
done  by  examining  IAW(I,  J),  If  there  is  an  RV  in  IAW  (IMIS, 
J)  the  subroutine  MCDWD  (Midcourse  Defense  with  Decoys)  is 
called,  with  calling  sequence  MCDWD  (IND, IAWT, IMIS,  J, 
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MNWT,  PMD,  PFA1,  PFA2,  PK,  INDSHT).  The  subroutine 
decreases  IAWT  (IMIS,  J)  if  an  RV  is  killed  and  sets 

INDSHT=0  when  no  midcourse  defender  expended 

INDSHT=1  when  midcourse  defender  expended  (either 
against  real  target  or  false  target). 

The  first  indicator  in  the  calling  sequence,  IND,  is  set  to  0 
for  random  defense  and  1  for  efficient  defense. 

The  subroutine  MCDWD  operates  as  follows.  The  total 
number  of  objects  NOB  JT  on  a  missile  is  counted,  made  up  of 
the  components: 

NOBJ1  =  number  of  RVs 

NOB  J2  =  number  of  heavy  decoys 

NOBJ2  =  number  of  light  decoys 

where 

NOBJT  =  NOBJ1  +  NOBJ2  +  NOBJ3. 

Each  RV  is  treated  as  though  it  comes  from  a  cluster 
characterized,  by  the  initial  mix  of  objects  as  above,  except 
when: 

IND=1,  NOBJ1  is  computed  from  IAWT. 

IND-0,  NOBJ1  is  computed  from  LAW 

The  group  of  objects  totaling  NOBJT  is  processed  as 
follows.  A  random  number  is  drawn  and  an  object  chosen 
which  matches  it.  (For  instance,  if  there  are  4  RVs,  2  heavy 
decoys  and  4  light  decoys  and  the  random  number  -  .05  the 
object  is  an  RV,  while  if  the  random  number  -  .91  the  object  is 
a  light  decoy. 

Another  random  number  A  is  chosen.  If  the  object  is  an  RV 
and  if  A  <  PMD  (probability  of  missed  detection)  it  is  not  shot 
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at  A  defender  is  not  wasted  but  neither  is  an  RV  engaged.  If 
A  >  PMD  the  RV  is  engaged  by  drawing  a  random  number  A 
and  if  A  £  PK  the  RV  is  killed 

If  the  object  is  a  decoy  PFA  is  set  equal  to  PFA1  or  PFA2, 
depending  on  the  type.  A  random  number  is  chosen.  If  A  <, 
PFA  there  is  a  false  alarm  and  midcourse  defender  is  wasted, 
e.  Decovs.  Efficient  Defense 


This  portion  of  the  program  is  analogous  to  the  no  decoys, 
efficient  defense  portion  discussed  in  b.  above.  However,  the 
subroutine  MCDWD  is  called. 

Before  entering  MCDWD,  INI>=1,  so  the  number  of  RVs  is 
computed  from  IAWT  rather  than  LAW  as  in  the  above 
discussion.  Otherwise,  the  procedure  is  the  same  as  discussed 
in  d.  above. 


f.  Decoys.  Preferential  Defense 

Not  currently  coded. 


7.  Subroutine  TD  fNVT.  NMIS.  ^ 


1EIL.IND) 


Subroutine  TD  (NVT,  NMIS,  MNWT,  PK IND),  Terminal 
Defense,  receives  from  MAIN  the  number  of  value  targets  NVT, 
number  of  missiles  in  the  attack  NMIS,  maximum  number  of  value 
targets  attacked  by  any  missile  MNWT,  probability  of  kill  PK,  and 
indicator  of  defense  option  IND. 

The  number  of  warheads  arri  ving  at  each  target  is  calculated  for 
all  targets,  namely  NW  (MAIN),  MAIN  =  1,  NVT.  Then  for  all 
targets  one  of  four  attrition  methods  is  utilized  in  determining  the 
number  of  warheads  surviving  after  the  terminal  defense,  or  the 
final  status  of  NW  (IMA IN),  MAIN  =  1,  NVT. 
a  Preallocated  Fixed-Salvo  Terminal  Defense 


A  particular  doctrine  is  chosen  for  the  targets  defended  by 
each  number  of  interceptors.  In  tire  example,  there  are  options 


for  2, 5  and  10  interceptors  at  each  target.  Their  firing 
doctrines  are,  respectively, 

2  interceptors  —  1,1 
5  interceptors  —  2,  1,  1, 1 
10  interceptors  --  2,  2,  2,  1,  1,  1,  1 
For  instance,  with  five  interceptors,  2  are  fired  at  the  first 
RV,  1  at  the  second,  1  at  the  third,  1  at  the  fourth,  and  none 
thereafter.  This  defense  should  be  conceptualized  as  fitting  a 
situation  where  the  attack  is  sequential  and  of  size  unknown  to 
the  defender.  For  10  interceptors,  the  "price’'  of  the  target  is  8, 
which  surely  kills  it  The  defense  is  balanced  so  that  the 
attacker  has  about  the  same  expected  value  per  RV  up  to  an 
attack  size  of  8  RVs.  This  preallocated  fixed-salvo  defense  is 
analogous  to  the  "Prim-Read"  defense. 

For  each  target  IMA  IN  =  I,  NVT,  the  number  of  warheads 
NWT  destined  for  it  is  determined.  The  number  of 
engagements  JENGS  is  chosen  corresponding  to  the  firing 
doctrine  associated  with  the  number  of  interceptors  defending 
the  target.  For  RV  1  through  NWT  the  number  of  interceptors 
fired  is  chosen  from  the  firing  doctrine.  As  RVs  are  killed  they 
are  subtracted  from  NW  (IMAIN). 

b.  Efficient  Terminal  Defense 

At  each  target  indexed  by  IMAIN  =  1,  NVT,  there  are  NW 
(IMAIN)  RVs  arriving.  There  are  1TD  (IMAIN)  terminal 
interceptors.  Set  NWT  =  NW  (IMAIN)  and  ITDT  =  ITD 
(IMAIN). 

If  ITDT  <,  NWT,  there  are  equal  to  or  fewer  interceptors  than 
RVs.  A  one-on-one  engagement  is  simulated  against  ITDT  of 
the  RVs.  The  number  of  successful  RVs  is  noted. 

If  ITDT  >  NWT  there  are  more  interceptors  than  RVs.  As  in 
the  efficient  boost-phase  defense  and  efficient  midcourse 
defense  previously  discussed,  there  are  IRP1.  defenders  shot  at 


.  a  fraction  PI  of  the  RVs  and  IR  defenders  shot  at  a  fraction  (1- 
Pl)  of  the  RVs. 

Each  RV  killed  causes  NW  (TMAIN)  to  be  decremented  by 
one.  The  final  number  of  successful  RVs  is  the  remainder  in 
NW  (1MA1N)  after  the  terminal  defense. 

c.  Limited  Shoot-Look-Shoot  Terminal  Defense 

In  this  case  the  defender  need  not  fire  salvos  at  RVs.  He  can 
shoot  an  interceptor,  observe  if  the  RV  is  killed,  shoot  again  if 
necessary,  and  so  on,  up  to  a  certain  limit  per  incoming  RV. 
Thus  he  can  substantially  conserve  defenders,  and,  if  the 
upper  limit  on  number  of  shots  is  reasonably  high,  have  a  high 
probability  of  kill  against  each  RV. 

The  program  is  se  t  to  have  the  upper  limit  of  the  truncation 
limit  of  the  minimum  of  Value/5  and  2.  In  the  example,  the 
value  of  the  attacked  target  is  either  20  or  10,  so  there  are 
effectively  4  or  2  shots  as  an  upper  limit 

The  number  of  terminal  defenders  is  ITDT.  The  number  of 
warheads  heading  for  the  target  is  NWT.  The  maximum 
number  of  shots  is  MAXDEF.  For  RVs  1  through  NWT, 
interceptors  are  fired  one  by  one  up  to  a  maximum  per  target  of 
MAXDEF.  If  an  interceptor  is  successful,  as  determined  by 
the  random  number  drawn  being  equal  to  or  less  than  PK,  then 
NW (TMAIN)  is  decremented  by  one.  As  each  interceptor  is 
fn.ied  ITDT  is  decremented  by  one.  The  process  terminates 
when  no  more  RVs  are  coming  or  when  ITDT  is  0.  The  final 
number  of  successful  RVs  is  the  remainder  in  NW(IMAIN). 

d.  Unlimited  Shoot  Look-Shoot  Terminal  Defense 

The  logic  for  this  defense  is  exactly  the  same  as  c.  above, 
except  there  is  no  upper  limit  on  the  number  of  shots  taken  at 
any  particular  RV.  There  is  no  check  of  MAXDEF  as  in  the 
above  discussion. 
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e.  Relative  Frequency. Array,  of  Successful  Penetrators 

At  the  end  of  Subroutine  TD  a  relative  frequency  array  of 
successful  penetrators  is  created.  The  array  is  denoted 
NWRF. 

In  the  example  there  may  be  any  number  from  0  through  20 
successful  penetrators,  so  the  relative  frequency  array  NWRF 
has  21  elements 

8.  Subroutine  VALSURV  (NVT.  IVSS) 

Subroutine  VALSURV  (NVT,  IVSS)  receives  from  MAIN  the 
number  of  value  target  NVT  and  returns  the  value  sum  surviving 
IVSS.  It  simply  sums  the  values  IV(I),  I  =  1,  NVT  associated 
with  the  targets  receiving  zero  successful  penetrators  (those  with 
NW  (I)  =0, 1=1,  NVT). 


EXAMPLE 


This  example  is  presented  in  order  to  give  the  .eader  information  on  the  level  of 
detail  at  which  the  model  is  used  and  of  the  sensitivities  of  the  outcome  to  variations  in 
BM/C^  options  and  defensive  resources.  Most  of  the  features  of  the  program  are  exercised 
and  illustrated. 

1.  Set  up  Targets  of  Differing  Value 

The  target  data  base  is  composed  of  the  array  IV  of  targets  of  differing 
value.  As  an  example,  consider  the  following  list  of  targets.  There  are  eleven 
conceptual  types  of  targets,  with  value  per  target  as  shown  and  total  value  as 
shown. 


Type  . of  Tarsal 

Number  of  Targets 

Value 

Per  Target 

TotaLYalys 

SLBM  Base 

3 

200 

600 

Installation 

10 

100 

1000 

ABM  Base 

10 

50 

500 

Bomber  Base 

50 

20 

1000 

Air  Defense  Base 

20 

20 

400 

CONUS  Military  Inst 

200 

20 

4000 

Theater  Military  Inst. 

200 

20 

4000 

Installation 

40 

10 

400 

CONUS  Econ.  Inst 

500 

10 

5000 

Theater  Econ.  Inst 

500 

10 

5000 

ICBM 

1QQQ 

5 

5Q QQ 

2,533 

26,900 

The  target  data  base  illustrated  consists  of  2533  targets  with  a  total  value  of 
26,900. 

The  subroutine  VALUES  creates  the  target  data  base. 

2.  Set  Up  Terminal  Interceptors 

For  each  target,  a  certain  number  of  terminal  interceptors  is  deployed  in 
array  ITD.  As  an  example,  for  each  target  in  the  data  base  given  abo^e, 


interceptors  per  target  are  deployed  equal  to  half  of  the  value  per  target  (rounded 
down  to  the  next  lower  integer).  The  result  is  as  follows: 


Type  of  Target 

Number 
of  Targets 

Value 

Per  Target 

Interceptors 

Eer  Target 

Total 

Interceptors 

5LBM  Base 

3 

200 

100 

300 

Installation 

10 

100 

50 

500 

ABM  Bases 

10 

50 

25 

250 

Bomber  Bases 

50 

20 

10 

500 

Air  Defense  Bases 

20 

20 

10 

200 

CONUS  Military  Inst. 

200 

20 

10 

200 

Theater  Military  Inst 

200 

20 

10 

2000 

^Installations 

40 

10 

5 

200 

CONUS  Econ.  Inst 

500 

10 

5 

2500 

Theater  Econ.  Inst. 

500 

10 

5 

2500 

ICBMs 

ioqq 

5 

2 

2000 

2,533 

12,950 

There  arc  12,950  terminal  interceptors  assigned  to  the  2533  targets.  The 
subroutine  TDINV  creates  the  terminal  interceptors  as  signments  in  array  ITD. 

3.  Set  Up  Attack  Allocation 

The  example  attack  allocation  consists  of  attacking  the  first  500  targets 
defended  by  fewer  than  25  interceptors,  namely  targets  24  through  523.  Each 
target  is  attacked  by  the  RVs  from  three  missiles.  Target  24,  for  instance,  is 
attacked  by  missile  number  1  with  5  RVs,  by  missile  501  with  5  RVs,  and  by 
missile  number  1001  with  10  RVs.  The  attack  is  as  follows: 

Missile  Number  Target  Number  K.Y2.££r_Iargei 

i  24  5 

500  523  5 

501  24  5 
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1000 

523 

5 

1001 

24 

10 

1500 

523 

10 

The  above  example  includes  cross-targeting  in  the  assignment  of  more 
than  one  missile  to  each  target  The  computer  program,  however,  is  more 
flexible  in  the  sense  that  the  logic  permits  each  RV  to  be  cross-targeted.  For 
instance,  missile  1  could  logically  fire  1  RV  each  at  targets  24, 25, 26, 27  and 
28.  (The  program  dimensions  are  currently  limited  to  the  targeting  at  most  two 
targets  by  each  missile,  but  the  logic  is  completely  flexible). 

From  the  point  of  view  of  the  targets,  the  above  allocation  results  in  the 
following  assignment  of  RVs: 


Terminal 

Attacking 

Target 

Value 

Interceptors 

RVs 

Number 

EfiLlargSt 

Per  Target 

Per  Target 

24 

20 

10 

20 

493 

20 

10 

20 

494 

10 

5 

20 

523 

10 

5 

20 

The  salient  point  is  that  each  target  receives  20  RVs,  from  three 
missiles.  Since  the  boost-phase  defense  kills  missiles  before  they  release  RVs 
each  target  can  be  attacked  by  0,  5, 10, 15  or  20  RVs,  depending  on  which 
missiles  (if  any)  ate  destroyed. 

Another  interpretation  worth  keeping  in  mind  is  that  the  example  attack 
allocation  might  be  conceptualized  as  assigning  1000  ICBMs  with  5  RVs  each, 
plus  500  SLBMs  with  10  RVs  each,  to  5C0  relatively  high  value  targets. 
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4.  R&&lis_foi:  12^Q.Xerminftl  Interceptors 

Figure  1  presents,  on  the  left  side,  results  for  1000  boost- phase 
interceptors  and,  on  the  right  side,  results  for  2000  boost-phase  interceptors. 
Units  of  value  destroyed  are  presented  as  a  function  of  number  of  midcourse 
interceptors  for  a  wide  variety  of  BMC-’  options. 

On  tiie  left  side,  graph  A, A, A  is  foi  random  allocation  of  boost-phase, 
random  allocation  of  midcourse,  and  preallocated  fixed-salvo  defense  at  each 
target.  It  yields  the  most  damage  for  the  attack  of  10,000  RVs.  Graph 
A,A,D,  with  unlimited  shoot-look-shoot,  does  quite  a  bit  better,  and  damage  is 
relatively  small  with  5000  midcourse  interceptors. 

Graphs  B,B,A  and  B,B,D  display  much  more  defense  effectiveness  than 
A, A, A  and  A,A,D.  In  particular,  B,B,D  with  3000  interceptors  yields 
relatively  small  damage.  Note  that  B,B,D  as  compared  with  B,B, A,  and 
A,AJD  as  compared  with  A, A, A  are  significantly  better.  The  inventory  of 
12,950  terminal  defenders  with  unlimited  shoot-look-shoot  is  much  more 
effective  than  when  these  same  terminal  defenders  are  limited  to  a  preallocated 
fixed-salvo  firing  doctrine. 

Finally,  on  the  left  side,  consider  B,CJD.  It  is  the  best  defense  with 
1000  interceptors  but  does  not  do  much  better  than  B,B,D  with  3000 
interceptors.  (There  will  be  more  discussion  of  preferential  midcourse  defense 
in  the  following  section). 

Turning  to  the  right  side  of  Figure  1,  note  that  2000  efficient  boost- 
phase  interceptors  (B,B,A  and  B,B,D)  achieve  near- zero  damage  with  as  few 
as  1000  efficient  midcourse  interceptors.  Confronting  1500  missiles  with 
2000  boost  phase  interceptor  shots  efficiently  allocated,  and  confronting  the 
remaining  RVs  with  1000  midcourse  interceptors  efficiently  allocated,  results 
in  expected  damage  of  zero  for  B,B,D  and  a  very  small  amount  for  3,B,A. 

Staying  with  the  right  side  of  Figure  1,  consider  the  random  defenses 
A,A,A  and  A,A,D.  The  former  requires  5000  midcourse  interceptors  to 
achieve  near- zero  damage,  while  the  laiter  achieves  this  with  3000  nvdeourse 
interceptors. 
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Comparing  the  left  with  the  right  sides  of  Figure  1,  note  that  1000  boost-phase 
interceptors  and  3000  midcourse  interceptors  with  B,B,D  does  almost  as  well 
as  2000  boost-phase  interceptors  and  5000  midcourse  interceptors  with  A, A, A. 
The  latter  needs  twice  the  number  of  boost-phase  interceptors  and  almost  twice 
the  number  of  midcourse  interceptors  to  achieve  the  same  effectiveness.  This 
illustrates  the  importance  of  BM/C^  capabilities  for  weapon  allocation. 


5.  Results  for  12.950  and  5.950  Terminal  Interceptors 

The  next  case  is  to  let  the  number  of  terminal  interceptors  per  target  be  equal 
to  one-fourth  of  the  value  per  target  (rounded  down  to  the  next  lower  integer). 
The  result  is  as  follows: 


T'vne  of  Target 

Number 
Per  Target 

Value 

Per  Target 

Interceptors 
Per  Target 

Total 

Intercep 

SLBM  Base 

3 

200 

50 

150 

C’3  Installation 

10 

100 

25 

250 

ABM  Base 

10 

50 

12 

120 

Bomber  Base 

50 

20 

5 

250 

Air  Defense  Base 

20 

20 

5 

100 

CONUS  Military  Inst  200 

20 

5 

1000 

Theater  Military  Inst. 

200 

20 

5 

1000 

C^  Installation 

40 

10 

2 

80 

CONUS  Econ.  Inst. 

500 

10 

2 

1000 

Theater  Econ.  Inst. 

500 

10 

2 

1000 

ICBMs 

1000 

5 

1 

1QQQ 

2,533 

5,950 

There  are  thus  5,950  as  opposed  to  12,950  terminal  interceptors. 

Figure  2  presents  the  results.  First,  on  the  left  side,  consider  the  solid  lines 
for  A,A,A  and  A,A,D.  The  two  cases  are  quite  close  together,  whereas  in  Figure  1 
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they  had  been  far  apart  (as  they  are  in  the  dashed  lines  of  Figure  2).  When  there 
are  fewer  terminal  interceptors,  whether  they  are  preallocated  or  shoot-look-shoot 
does  not  yield  so  much  difference  in  payoff. 

Second,  on  the  left  side,  consider  the  solid  lines  for  B,B,A  and  B,B,D. 
There  is  a  dramatic  impact  between  3000  and  4000  midcourse  interceptors.  The 

payoff  from  taking  on  all  of  the  RVs  surviving  the  boost-phase  defense  is  immense. 
When  there  are  1000  boost-phase  interceptors  with  of  0,9,  there  are  roughly  600 

surviving  missiles  of  which  400  have  5  RVs  and  200  have  10  RVs,  or  4000  RVs. 

If  these  can  all  oe  engaged  the  number  of  survivors  is  roughly  400,  distributed 
(quite  non-uniformly)  across  the  500  attacked  targets.  Both  types  of  terminal 
defenses  are  effective,  with  the  un’: united  shoot-look-shoot  doing  better.  The  main 
point,  however,  is  that  having  enough  efficient  midcourse  defenders  to  engage  all 
of  the  RVs  makes  a  tremendous  difference. 

Third,  on  the  left  side,  consider  the  solid  line  for  B,C,D.  Up  to  3000 
midcour&e  defenders  preferential  midcourse  defense  does  quite  a  bit  better  than 
efficient  midcourse  defense  (B,B  J)).  Beyond  this  the  algorithm  does  not  beat  the 
uniform  defense.  If  the  defender  were  to  observe  the  attack  size,  he  could  simply 
switch  to  uniform  defense  if  he  knew  the  attack  size  exceeded  4000. 

In  general,  the  preferential  defense  B,C,D  is  a  difficult  optimization 
problem.  The  algorithm  described  previously  is  independent  of  attack  size.  A 
better  algorithm  would  solve  the  large-  scale  integer  assignment  of  shooters  to 
targets  knowing  both  numbers  of  shooters  and  targets,  and  knowing  the  terminal 
defense  capabilities.  However,  the  present  example  makes  the  point  that 
preferential  midcourse  defense  is  much  better  than  efficient  rrudeourse  defense 
when  the  number  of  midcourse  defenders  is  limited. 

Fourth,  consider  the  right  side  of  Figure  2.  As  before,  if  there  are  2000 
efficient  boost -phase  interceptors  and  1000  or  mere  efficient  midcourse  interceptors 
(B,B,A  or  B,B,D)  there  is  very  little  damage  for  the  solid  or  dashed  lines.  Halving 
the  numbers  of  terminal  interceptors  is  not  a  problem.  However,  if  both  the  boost- 
phase  and  midcourse  interceptors  are  random  (A, A, A  and  A,A,D)  the  solid  lines  are 
significantly  worse  than  the  dashed  ones. 


In  particular,  A,A,D  is  much  worse,  for  unlimited  shoot,-look-shoot  with 
fewer  terminal  interceptors  is  negated  by  exhaustion  at  quite  a  few  targets.  Where 
3000  midcourse  interceptors  result  in  near- zero  damage  with  the  dashed  line, 

A, A  JD,  5000  midcourse  interceptors  are  required  with  the  dashed  line  A,A,D. 

EfTicient  Bocst-Phase,  Random  Midcourse 

The  right  side  of  Figure  3  presents  results  for  efficient  boost-phase  defense 
followed  by  random  midcourse  defense.  The  dotted  lines  are  the  new  results. 

There  are  few  enough  missiles  surviving  the  efficient  boost-phase  defense 
that  random  midcourse  defense  does  not  result  in  significantly  more  damage  than 
efficient  midcourse  defense. 


Random  Boost-Phase.  Efficient  Midcourse 

The  right  side  of  Figure  4  presents  results  for  random  boost-phase  followed 
by  efficient  midcourse  defense.  The  dotted  lines  are  the  new  results. 

Note  that  2000  random  boost-phase  defenders  let  through  so  many  RVs  that 
1000  efficient  midcourse  interceptors  do  not  do  much  good. 

However,  as  the  numbers  of  efficient  midcourse  interceptors  increases  from 
2000  to  3000  the  attack  is  sufficiently  diluted  so  that  just  a  few  targets  are 
destroyed.  This  midcourie  defense  situation  is  exactly  analogous  to  the  left  side  of 
the  figure  increasing  from  3000  to  4000  efficient  midcourse  interceptors. 

■  Random  Bosst-Phasc.  Hybrid  Mdcourec 

The  right  side  of  Figure  5  presents  results  for  random  boost-phase  defense 
followed  by  a  combination  of  efficient  midcourse  defense  and  random  midcourse 
defense.  The  dotted  lines  are  the  results. 

Hybrid  defense  is  defined  as  one-half  efficient  defense  followed  by  one-half 
random  defense.  (It  is  denoted  by  B'  =  1/2B,  1/2A).  Thus,  the  same  number  of 
RVs  after  boost-phase  defense  as  in  Figure  4  is  confronted  by  a  partially  efficient 
and  partially  random  midcourse  defense. 


Figure  3.  EFFICIENT  BOOST-PHASE,  RANDOM  MIDCOURSE 
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Figure  4.  RANDOM  BOOST-PHASE,  EFFICIENT  MIDCOURSE 


Figure  5.  RANDOM  BOOST-PHASE,  HYBRID  MIDCOURSE 


The  general  behavior  is  for  the  units  of  value  destroyed  to  decrease  more 
gradually  as  a  function  of  number  of  midcourse  defenders.  The  number  of 
midcourse  defenders  to  yield  low  damage  increases  from  3000  with  efficient 
defense  (A,B,A  and  A3JD)  to  4000  with  hybrid  defense  (A,B',A  and  A,B',D). 

9.  Hybrid  Boost-Phase.  Hybrid  Midcourse 

The  right  side  of  Figure  6  presents  results  for  hybrid  boost-phase  defense 
followed  by  hybrid  midccurse  defense. 

In  this  case,  a  hybrid  boost-phase  defense  is  composed  of  1000  efficient 
boost-phase  defenders  and  1000  random  boost-phase  defenders.  The  1000 
efficient  defenders  engage  a  subset  of  the  1500  missiles  and  the  1000  random 
defenders  engage  the  remaining  missiles  (roughly  600  missiles).  Thus  the  number 
of  RVs  proceeding  to  be  confronted  by  the  midcourse  defense  is  more  than  in  the 
efficient  boost-phase  case  but  less  than  in  the  random  boost-phase  case. 

As  expected,  the  effects  of  adding  hybrid  midcourse  interceptors  are 
somewhat  gradual,  with  value  destroyed  being  fairly  low  with  3000  defenders. 

The  gradual  effects  of  the  hybrid  midcourse  defense  can  be  comrasted  with  the 
sharp  effects  of  the  efficient  midcourse  defense  (see,  for  example,  B',B',A  and 
B'JB'JD  versus  A,B,A  and  A,B,D). 

10.  Hybrid  Boost-Phase.  Random  Midcourse 

The  right  side  of  Figure  7  presents  the  last  excursion  of  this  type. 

Hybrid  boost-phase  defense  is  followed  by  random  midcourse  defense. 
The  effects  of  adding  random  midcourse  defenders  are  more  gradual  than  adding 
hybrid  midcourse  defenders. 

11.  Suppression  of  Boost-Phase  Defense 

The  right  side  of  Figure  8  presents  results  for  two  types  of  suppression  of 
boost-phase  defenses. 

In  the  first  type  of  suppression  boost-phase  defenders  are  reduced  by  one- 
third,  or  from  2000  to  1333.  Not  all  of  the  missiles  are  engaged.  When  the 
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Figure  6.  HYBRID  BOOST-PHASE,  HYBRID  MIDCOURSE 


Figure  7.  HYBRID  BOOST-PHASE,  RANDOM  MIDCOURSE 
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Figure  8.  SUPPRESSION  OF  BOOST-PHASE  DEFENSE 


number  of  efficient  midcounie  defenders  is  sufficient  to  engage  ail  of  the  surviving 
RVs,  units  of  value  destroyed  is  reduced  to  a  small  amount.  Again,  the  behavior  of 
results  is  the  same  as  in  previous  examples  with  sufficient  numbers  of  midcourse 
defenders. 

In  the  second  type  of  suppression  one- third  of  the  missiles,  at  random,  are 
not  engaged.  The  number  of  efficient  midcourse  defenders  required  to  result  in 
small  damage  is  larger  than  in  the  first  case. 

12.  Preferential  Targeting  by  Boost-Phase  Defenders 

The  left  side  of  Figure  9  presents  results  for  boost-phase  defensive  firing 
proportional  to  number  of  RVs  on  the  missiles  being  attacked. 

Two  cases,  C,B,A  and  C.BJD,  are  presented.  They  can  be  compared  with 
B,B,A  arid  B.BJD,  and  are  seen  to  have  less  damage. 

The  reason  why  the  return  to  the  defender  is  not  greater  for  proportional 
targeting  of  boost-phase  defenders  is  that  the  attacking  missiles  have  either  5  or  10 
RVs  each.  If  tthe  mix  of  RVs  were  from  1  to  20  per  missile,  or  something  quite 
skewed,  the  improvement  for  the  defender  would  be  much  greater. 


The  left  side  of  Figure  10  shows  the  effects  of  three  decoy  situations. 


To  begin,  replace  two  of  each  five  RVs  by  four  heavy  decoys,  resulting  in  a 
total  attack  consisting  cf  6000  RVs  and  8000  decoys.  In  all  cases,  assume  that  the 
BM/C^  system  of  the  defense  is  B,B,A. 

In  the  first  example,  K  =  2.5,  which  can  be  interpreted  as  meaning  that  the 
inherent  error  of.  the  measuring  device  is  such  that  RVs  and  decoys  yield 
observations  which  are  two  and  one-half  standard  deviations  apart  A  detection 
threshold  is  set  resulting  in  PMD  =  .  1  and  PFA  =  .  1  (probability  of  missed 
detection  against  an  RV  equal .  1  and  probability  of  false  alarm  against  a  decoy 
equal  .1).  In  this  case,  the  resulting  units  of  value  destroyed  as  a  function  of 
number  of  midcourse  interceptors  is  lower  than  the  baseline  curve  B,B,A.  So  this 
is  not  a  good  attacker  option. 


Figure  9:  PROPORTIONAL  TARGETING  BY  BOOST-PHASE  DEFENSE 


In  the  second  example,  K  =  1.0,  and  the  detection  threshold  yields  PMD  = 
.3  and  PFA  =  .3.  This  is  much  worse  for  the  defense  than  K  =  2.5.  But,  still,  it 
is  not  a  good  attacker  option. 

In  the  third  example,  however,  assume  that  30  light  decoys  can  be  added 
with  PFA  =  .  1  (smaller  than  the  PFA  =  .3  for  heavy  decoys).  The  damage  is 
smaller  than  B,B,A  until  the  number  of  midcourse  interceptors  is  greater  than 
3000,  but  then  it  is  not  reduced  to  near- zero.  Even  with  5000  interceptors,  there  is 
significant  damage. 

The  third  example  is  a  striking  demonstration  of  the  difficulty  ir  to  the 
defense  of  dealing  with  decoys.  Near- zero  damage  becomes  very  difficult  to 
attain. 


14.  Concentrated  Attack 

This  example  is  very  different  from  those  previously  discussed,  for  in  all  of 
those  cases  the  attack  was  held  constant.  Now,  the  effects  of  concentrating  on  a 
smaller  set  of  targets  are  investigated. 

Allocate  all  of  the  missiles  to  the  50  targets  24  through  73.  The  maximum 
possible  amount  of  value  destroyed  is  1000  units  as  compared  to  9700  units 
previously. 

Both  the  left  and  right  sides  of  Figure  1 1  show  results  of  the  new  attack. 

On  the  left  side  the  BM/C^  assumption  is  B,B,D  while  on  the  right  side  it  is 
A,B,D. 

On  the  left  side,  for  the  original  attack,  4000  interceptors  yields  near- zero 
damage  while  for  the  near  attack  even  6000  interceptors  result  in  significant 
damage.  'Riis  is  because  there  are  200  RVs  initially  aimed  at  each  target,  and  the 
process  is  sufficiently  "lumpy"  that  at  a  few  targets  there  are  enough  RVs 
remaining  to  exhaust  the  terminal  defense. 

On  the  right  side,  for  the  original  attack  3000  efficient  interceptors  were 
enough  to  significantly  blunt  the  attack.  Here  4000  are  required. 
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Figure  11.  CONCENTRATED  ATTACK 


It  is  desirable  to  compare  the  results  of  "expected  value"  calculations  with 
Monte  Carlo  calculations  in  oroer  to  illustrate  the  large  differences  which  can 
occur  under  certain  circumstances.  This  is  done  below  for  two  cases. 

In  both  cases  the  numbers  of  defenders  are  as  follows: 

Boost-Phase  1500 
Midcourse  2000 

Terminal  5950 

The  expected  value  calculation  procedure  follows.  Assume  that  1500 
efficiently  allocated  boost-phase  defenders  confront  1500  missiles,  and  the  is 

.8.  The  expected  number  of  missiles  surviving  is  300,  comprised  of  200  missiles 
with  5  RVs  and  100  missiles  with  10  RVs,  for  a  total  of  2000  RVs.  Assume  that 
2000  efficiently  allocated  midcourse  defenders  confront  the  2000  RVs.  The 
expected  number  of  survivors  is  400. 

Now  assume  that  400  RVs  are  uniformly  distributed  over  the  500  targets. 
Since  there  are  470  targets  of  value  20  and  30  targets  of  value  10  assume  that  376 
of  the  first  type  and  24  of  the  second  type  are  attacked  by  one  R  V  each.  With 
preallocated  fixed-salvo  defease  of  2  interceptors,  the  probability  of  penetration  is 
.04,  so  376  x  20  x  .04  =  300.8  units  of  value  are  killed  at  the  higher-value 
targets.  Similarly,  24  x  10  x  .2  =  48  units  of  value  are  killed  at  the  lower-value 
targets.  The  total  is  348.8.  With  efficiently  allocated  terminal  defense  the 
equivalent  figures  are  376  x  20  x  .00032  =  2.406  and  24  x  10  x  .04  -  9.6.  The 
total  is  12.006. 

Running  the  Monte  Carlo  model  for  30  trials  yields,  for  the  two  cases 
above,  ihe  foliowing  results: 

MaoJlLSanipte  Standard  Deviation  of  Sample 
Preallocated  Terminal  Defense  1184  286 

Efficient  Terminal  Defense  844  344 

The  inference  about  the  process  from  the  sample  is  as  follows: 


i 


Standard  Deviation  of 
Estimate  cf  Population 
Estimate  of  Popu  lation  Mean  Mean 

Preailocated  Terminal  Defense  1184  53 

Efficient  Terminal  Defense  844  64 


Summarizing  the  expected  Monte  Carlo  value  calculations  the  following  is 
obtained: 


Estimate  of 

Estimate  of 

Standard  Deviation  of 
Estimate  of 

Mean  with 

Mead  with 

Mean  with 

Expected  Value 

Monte  Carlo 

Monte  Carlo 

Cal'iaMon 

Calcuialmn 

Preallocated  Terminal  Defense 

349 

1184 

28 

Efficient  Terminal  Defense 

12 

844 

54 

Performing  exactly  the  same  analysis  with  probability  of  kill  equal  .9  yields 
the  following  (40  Monte  Carlo  trials): 


Estimate  of 
Mean  with 
Expected  Value 
Q&uk&m 


Estimate  of 
Mean  with 
Monte  Carlo 
Cakatfation 


Standard  Deviation  of 
Estimate  of 
Mean  with 
Monte  Carlo 
Cakuktion 


Preailocated  Terminal  Defense  25  246  28 

Efficient  Terminal  Defense  1  126  20 


Fignt/e  12  and  13  show  the  calculations  leading  to  tire  above  results. 


a 


I 


m 


% 


■m 

tl 


y>71 


OU 

CO 

30 

i! 

M 

II 

jtf 

J3C 

a. 

CL 

Q. 

2 

2 

2 

H" 

{— 

BS 

3B 

35 

«o 

60 

OO 

CC 

oe 

CC 

LU 

1X1 

UJ 

© 

© 

© 

wr 

U&. 

2 

2 

UJ 

Ui 

UJ 

u_ 

u„ 

i » 

UI 

UJ 

UJ 

© 

© 

© 

© 

o 

© 

o 

© 

m 

in 

© 

cn 

T— 

CM 

in 

uu 

^  w 

uu  ~ 

u,  jf- 

ui  UJ 

■“  ijj 
ULi  O 

<  tu 
■E  9/9 

a.  SE 

i  j 

f—  o 
eo  o 
o  o 
o  = 

03  2fc 


X 

o 

o 

in 

. .  + 

uj  m 


+ 

>n 

X 

o 

O  . 
CM  w 
^  O 

ui  ^ 
^  ui 

iu 

.  i  *— 

UJ  Ui 

o  }h 

M  {jj) 

Mi  © 

5  CO 

4-  cr 
gv.  “ 

.  © 

5—  O 

60  « 
©  © 
o  ««• 
03  £ 

os  c c 

UI  UJ 
h"  I— 
u.  Lu 
<  < 

CO  CO 
>  2> 
re  ce 

u_  uu 

o  © 

HC  CC 
UJ  Ui 

ca  co 
£  £ 


co 

4. 

UU  ^ 

©  r~, 

CC  CO 

<e 

t  o 

o 

©  CO 

©  I! 

os*  " 

©  © 
T  X 


®  (gy 

©  ?r 
^  w 
uu  ~“i 
°  © 

2  * 

©  H 

r~  otf 

&  *• 

O  2 
©  J— 

d  5 

*  Ui 
“g  CO 
£  2r- 
O  £ 

”  UJ 
Z  Q 

=>  -j 

w  <5 

z  Sn 


,  UJ 
~  m  CU  6V 
<C  2  uj 

2  CO  ©  III 
£  ?l!  d  S’ 


»  us 

«  a  II 

u,  UJ  CM 
ULI  1h  * 

«  <  x 

£  o  T- 


<  « 
35  © 

£  ,! 
©  S 

ui  °. 

H-  x 

►”  © 
SB  r- 


X 

s 

u-  ac 

o  ** 

*** 

*3**  *•» 

©Sir 

iSi 

25  Si 

OO  b 
«  *-  !l 

Cg  LU  . — . 

<  ®  3 

®  Ui  to 

2w>  kw> 

<  < 


co  'cr 
in  us 


•*£ 

sa 


Q  ‘i*  S£  2 

<7Q  ^ 

*  O  to  {V  CO 

co  < 
o 


II  - 

ix  o 
35  uj  uj 
£5%2e 

uj  k-  S2 

jl  W  H 

su5 

|oa 

w  <  o 
a. 


C 

O  h- 
O  2 


^  O 

2  X 
2  © 
£  c 
ec  !— • 
uj  *— 
h—  «s£ 


© 

© 

< 

'HI  X 

UJ 

X 

CO 

c 

— j 

UJ 

2 

2 

© 

•i  . 
<  ^ 

u 

a 

an] 

< 

U 

O 

£31 

© 

Ui  cu 

LU 

o 

© 

LU 

U. 

Ui 

UJ 

Ui 

CC 

U_ 

— 1 

LU 

CC 

LU 

h” 

►  - 

(— 

ex, 

LU 

H~ 

CU 

UJ 

© 

u 

o 

<c 

CJ 

Ui 

Ui 

Ui 

ui 

33k. 

a. 

cu 

u 

cu 

X 

X 

X 

LU 

X 

Ui 

Ui 

UJ 

kM 

H 

UJ 

SW'i’iVWiv 


V- 1 


•i'4 

'*  rit 

v.l’t 


U/tf 


;<#S| 

p 


•..rv 


..v. 

•v  >.i 


m 


Vftl 


1 


-,W 

il  >V- 

'.'•■it,,. 


I 


m 


m 


v«F»j 

V*! 

■* . 


m 

m 


& 


CO 

X  • 

O  -5  «,  o 

J.  ^  CM  CM 


- —  tS 


o  o 
uu 

f3 60  2  tr  ^  p- 

2K  o  's'  £! 

aSz 


00 

t3 

%  + 

«3!  p~i 

H*  ee 


o  Ui 


«tt*jLL 


o  ®o 

O  ▼“ 

*”  I! 


O'  =1 

w2_ 
«£  U 

IU  Z 


CO  CO 
^  CM 
CM  r- 


1=  P 


an 

a> 

an 

li 

II 

II 

.XL 

4C 

•X 

a. 

CL. 

DL 

X 

X 

X 

1— 

>— 

H 

=S 

£ 

=£ 

60 

CO 

CO 

or: 

X 

ae 

UJ 

UJ 

LU 

0 

0 

a 

X 

X 

X 

UJ 

UJ 

LU 

LL. 

u. 

u. 

UJ 

UJ 

UJ 

O 

0 

0 

O 

0 

0 

O 

0 

in 

uo 

0 

cn 

T" 

m 

UJ 

5  (jj 

LL  *> 


4®  x 

«  s 

O  X 

w  m&> 
T— 


WJ  y 

Ct  H*i 


“#r  • 

o  » 


£•  st.  **  t>- 

s  -pr  cm 

«M,  O 

LCJ 

a 


“■  it* 

UJ  Q 

<  UJ 

2  0*3' 

£  SE 

I  =3 
O 
60  O 

O  O 

0  55 
oa  SE 


cj  x 


"?i< 


O  1. 


-*j  3s 
c 

UJ 

5E  wo 
ec  « 


u?  ca  -3-  co  jo 

S  x  .25,  cn  *»” 

X  <  O  ▼"  V" 


cc  as 

JU  UJ 


b»  C/J 

>  > 

CC  CC 


=?  UJ 

z  o 

‘■=>  „J 

u  ^ 

3?  X  »» 

X  “  CM 

S  1  III 


=i 

«  c  00  n 

Nt  Z  1^  t 
p  «£  CM  T— 


Q  CC 

2  lli 


w  "a 


X 

LU 

U. 

UJ 

a 

CO 

X 

LU 

LU 

LU 

UJ 

CO 

X 

LJ 

fU 

T 

in 

X 

0 

0 

w 

CC 

LU 

00 

SE 

w 

X 

UJ 

CD 

SE 

v. 

UJ 

CD 

< 

a 

UJ 

H 

<L 

O 

O 

LU 

O 

UJ 

T"* 

X 

_J 

CO 

c 

UJ 

CO 

cc 

0 

LU 

0 

X 

z 

0 

_J 

< 

JJ 

■mJ 

O 

Q 

0 

0 

UJ 

a. 

r* 

UJ 

LU 

UJ 

cc 

1 

X 

X 

X 

to 

H 

5— 

H 

a. 

H 

0 

0 

■*** 

O 

O 

O 

CO 

0 

SE 

cs: 

z 

LU 

LjJ 

LU 

0 

0 

t 

,IB— 

Cl. 

CU 

Q. 

0 

3! 

LU 

SE 

X 

X 

X 

C8 

c 

cc 

UJ 

UJ 

LJ 

cc 

II 

0 

MJ 

K* 

T" 

O 

•  • 

LU 

CD 

LU 

H 

c 

K 

X 

< 

O 

h* 

X 

<3 

£ 

O 

Z 

UJ 

r1 

< 

—J 

LU 

Cw 

X 

0 

— 1 
< 
UJ 

O 

LU 

CO 

0 

0 

LU 

LU 

_J 

LU 

cc 

LU 

LU 

at; 

t— 

CL 

WJ 

c 

0 

O 

UJ 

Cm 

UJ 

u. 

X 

UJ 

•  5  r  •  *  *  •-.  i  r  1 


i  i-innS  •  ■■ 


Figure  13.  C0MW^IS0Nr0F  “EXPECTED  VALUE'  AND  MONTE  CARLO 


DEFINITIONS  OF  VARIABLES 


Maine 

AVES 

AVED 

FIAWT 

FKLD 

FTARS 

IAD(I,J,K) 

lAT(IJ) 

IAW(I,J) 

IAWMT 

IAWS 

IAWT 

IAWT(U) 

IBPDS 


Definition 

Average  value  sum  surviving,  completed  cumulatively  as  of  each 
trial  Computed  and  used  in  MAIN. 

Average  value  sum  destroyed,  computed  cumulatively  as  of  each 
trial.  Computed  and  used  in  MAIN. 

Sum  of  total  number  of  RVs  in  array  IAW(I,J).  Computed  and 
used  in  Subroutine  BPD. 

Fraction  of  boost-phase  defenders  killed  in  boost-phase  defense 
suppression.  Set  by  user  at  beginning  of  Subroutine  BPD,  if 
desired. 

Fraction  of  missile  targets  which  cannot  be  attacked  by  boost- 
phase  defenders.  Set  by  user  at  beginning  of  Subroutine  BPD,  if 
desired. 

Number  of  decoys  of  type  K  in  the  Jth  group  of  RVs  carried  by 
missile  I.  In  COMMON. 

Target  aimed  at  by  all  of  the  RVs  in  the  Jth  group  of  RVs  carried 
by  missile  I.  In  COMMON. 

Number  of  RVs  in  the  Jth  group  of  RVs  carried  by  missile  I.  In 
COMMON. 

Sum  of  total  number  of  warheads  on  a  particular  missile. 
Computed  and  used  in  Subroutine  BPD. 

Sum  of  total  number  of  RVs  in  array  LAW(IJ).  Computed  and 
used  in  MAIN. 

Sum  of  total  number  of  RVs  in  array  IAW(I  J).  Computed  and 
used  in  Subroutine  BPD. 

Temporary  array  of  number  of  RVs  in  the  Jth  group  of  RVs 
carried  by  missile  I.  Computed  and  used  in  Subroutines  MCD, 
MCD3,  MCDWD. 

Indicator  for  boost-phase  defense  suppression.  If  IBPDS=0,  no 
boost-phase  suppression.  If  IBPDS=1,2,3  or  4, boost-phase 
suppression  of  different  types.  Set  by  user  at  beginning  of 
Subroutine  3PD. 


fflPDST 

IBPM 

ICMPLT 

DFRAC 

IMCM 

INDI(I,J) 

INDJ(I,J) 

INDTA(I) 

IRS 

itd® 

rros 

iv® 

IVS 

IV  SD 
IVSS 


Working  variable  indicator  for  boost-phase  defense  suppression. 
May  be  set  at  0  or  1.  Set  and  used  within  Subroutine  BPD. 

Indicator  for  hybrid  boost  phase  defense.  If  IBPM=0,  no 
hybrid.  If  IBPM=1,  hybrid.  Used  in  MAIN. 

Indicator  denoting  whether  boos  tiers  receive  complete  coverage 
(i.e.,  whether  all  boosters  necessarily  receive  at  least  one  shot  if 
possible).  i/ICMPLT=0,  no  complete  coverage.  If  ICMPLT=  1 , 
complete  coverage.  Computed  and  used  in  Subroutine  BPD. 

Indicator  for  fraction  of  defense  which  is  efficient  or  random. 
IFRAC-  T  denotes  efficient  portion  and  IFRAC=2  denotes 
random  portion.  Used  in  MAIN. 

Indicator  for  hybrid  midcourse  defense.  If  IMCM=0,  no  hybrid. 
If  IMCM=  1 ,  hybrid.  Used  in  MAIN. 

Indicator  for  the  Ith  target  and  the  Jth  missile  aimed  at  that  target, 
giving  the  number  of  the  missile.  Computed  and  used  in 
Subroutine  MCD3. 

Indicator  for  the  Ith  target  and  the  Jth  missile  aimed  at  that  target, 
giving  the  value  cf  J.  Computed  and  used  in  Subroutine  MCD3. 

Indicator  for  whether  or  not  Ith  missile  has  been  assigned  to  a 
boost-phase  defense  shooter.  If  INDTA(I)=0,  it  has  not.  If 
INDTA(I)=1,  it  has.  Computed  and  used  in  Subroutine  BPD. 

Seed  of  the  random  number  generator,  changed  each  time  a 
random  number  is  generated.  In  COMMON. 

Number  of  terminal  defenders  at  target  I.  In  COMMON. 

Sum  of  total  number  of  terminal  defends  s.  Computed  and  used 
in  MAIN. 

Value  of  target  I.  In  COMMON. 

Sum  of  total  number  of  units  of  value.  Computed  and  used  in 
MAIN. 

Value  sum  destroyed.  Computed  and  used  in  MAIN. 

Value  sum  surviving.  Computed  in  Subroutine  VALSURV. 

Used  in  MAIN. 


JPR(I) 

KPR(I,J) 

LPI 

NKLD 

NOBJ1 

NOBJ2 

NOBJ3 

NOBJT 

NTARS 

NVT 

NW(I) 

NWRF(I) 

NWRMP 

NWT(I) 

PFA1 


Number  of  shots  taken  by  preallocated  preferential  defense  with 
total  number  of  interceptors  I.  Specified  and  used  in  Subroutine 
TD. 

Number  of  interceptors  fired  at  the  Jth  incoming  RV  when  the 
total  number  of  interceptors  is  I.  in  a  preallocated  preferential 
defense.  Specified  and  used  in  Subroutine  TD. 

Long  printout  indicator.  If  LPI=0,  no  long  printout.  If  LPI=  J , 
long  printout  In  COMMON. 

Number  of  boost-phase  defenders  killed  in  boost-phase  defense 
suppression.  Set  by  user  at  beginning  of  Subroutine  BPD,  if 
desired. 

Number  of  objects  of  type  1,  namely  RVs.  Computed  and  used 
in  Subroutine  MCDWD. 

Number  of  objects  of  type  2,  namely  heavy  decoys.  Computed 
and  used  in  Subroutine  MCDWD. 

Number  of  objects  of  type  3,  namely  light  decoys.  Computed 
and  used  in  Subroutine  MCDWD. 

Number  of  objects  total.  Computed  and  used  in  Subroutine 
MCDWD. 

Number  of  missile  targets  which  cannot  be  attacked  by  boost- 
phase  defenders.  Set  by  user  at  beginning  of  Subroutine  BPD,  if 
desired. 

Number  of  value  targets.  Created  in  VALUES.  Passed  in  calling 
sequences  of  subroutines. 

Number  of  RVs  headed  for  target  I.  In  COMMON, 

Frequency  distribution  of  number  of  targets  receiving  I 
penetrators  at  the  end  of  the  simulation.  In  COMMON. 

Maximum  number  of  RVs  which  can  hit  a  target,  plus  1 .  In 
COMMON. 

Temporary  array  of  number  of  RVs  headed  for  target  I. 

Computed  and  used  in  Subroutine  MCD3. 

Probability  of  false  alarm  against  heavy  decoys,  where  midcourse 
defender  mistakenly  calls  a  decoy  an  RV.  Specified  in 
Subroutine  MCD  and  used  in  Subroutine  MCDWD. 


PFA2 

Probability  of  false  alarm  against  light  decoys,  where  midcourse 
defender  mistakenly  calls  a  decoy  an  RV.  Specified  in 

Subroutine  MCD  and  used  in  Subroutine  MCJDWO. 

PMD 

Probability  of  missed  detection,  where  midcourse  defender  does 
not  detect  an  RV.  Specified  in  Subroutine  MCD  and  used  in 
Subroutine  MCDWD. 

RAT 

Ratio  of  attacking  RVs  to  terminal  interceptors.  Computed  and 
used  in  Subroutine  MCD3. 

RATSL 

Lower  bound  on  ratio  of  attacking  RVc  to  terminal  interceptors 
for  particular  iteration  of  preferential  defense  algorithm. 
Computed  and  used  in  Subroutine  MCD3. 

RATSU 

Upper  bound  on  ratio  of  attacking  RVs  to  terminal  interceptors 
for  particular  iteration  of  preferential  defense,  algorithm. 
Computed  and  used  in  Subroutine  MCD3. 

SD 

Standard  deviation  of  value  sum  surviving,  computed 
cumulatively  as  of  each  trial.  Computed  and  used  in  MAIM. 

vssa) 

Working  variable  to  preserve  value  surviving  at  the  end  of  the  Ith 
trial  of  the  simulation.  Used  in  MAIN. 

APPENDIX  B 

LISTING  OF  COMPUTER  PROGRAM 
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